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ABSTRACT 

Research  Laboratory  Space  Time  Adaptive  Processing  (RLSTAP)  is  a  com¬ 
pilation  of  a  multiple-channel  radar  environment,  signal  and  hardware  pro¬ 
grams  integrated  within  the  KHOROS  environment.  RLSTAP  was  initially 
developed  for  the  US  DARPA  Mountain  Top  Program,  and  has  since  been  used 
by  several  TTCP  subgroups  for  collaborative  efforts.  Member  nations  of  the 
SEN  TP3  group  have  been  using  RLSTAP  as  the  baseline  radar  system  simu¬ 
lator  for  collaborative  radar  performance  assessment.  Under  SEN  TP3  Task  4 
Collaborative  Radar  Simulation  Validation,  we  have  integrated  an  L-band  sea 
clutter  model  with  existing  RLSTAP  programs  without  extensive  redesigning 
or  recoding  of  the  original  sea  clutter  “C”  program.  This  report  serves  to  illus¬ 
trate  how  the  KHOROS  software  development  tools  can  be  used  to  modify  a 
stand-alone  L-band  sea  clutter  model  written  in  “C”  for  KHOROS  compliance. 
The  rest  of  the  report  contains  an  overview  of  RLSTAP  clutter  models  and 
detail  formulation  of  the  second-order  L-band  sea  clutter  model. 


APPROVED  FOR  PUBLIC  RELEASE 


DEPARTMENT  OF  DEFENCE 

DEFENCE  SCIENCE  &  TECHNOLOGY  ORGANISATION 


(\Q 


DSTO-TR-1206 


Published  by 

DSTO  Electronics  and  Surveillance  Research  Laboratory 
PO  Box  1500 

Salisbury,  South  Australia,  Australia  5108 

Telephone:  (08)  8259  5555 
Facsimile:  (08)  8259  6567 

©  Commonwealth  of  Australia  2001 
AR  No.  AR-012-013 
September  2001 


APPROVED  FOR  PUBLIC  RELEASE 


DSTO-TR-1206 


Integrating  an  Airborne  L-band  Sea  Clutter  Model  into 
Research  Laboratory  Space  Time  Adaptive  Processing 

(RLSTAP) 


EXECUTIVE  SUMMARY 

This  work  has  been  performed  under  the  task  DST99/045- “RADAR  SYSTEMS  MOD¬ 
ELLING/PREDICTION  /VALIDATION”  for  TTCP  SEN-TP3  Task  4  Collaborative  Radar 
Simulation  Validation. 

Research  Laboratory  Space  Time  Adaptive  Processing  (RLSTAP)  is  a  compilation  of  a 
multiple-channel  radar  environment,  signal  and  hardware  programs  integrated  within  the 
KHOROS  environment.  The  KHOROS  software  development  tools  provide  a  complete 
environment  for  developing  new  toolboxes  and  programs.  These  programs  are  graphically 
represented  as  glyphs  in  the  KHOROS  environment  and  were  written  exclusively  for  mod¬ 
elling  monostatic  and  bistatic  radar  (RLSTAP  models  use  a  graphical  representation 
where  the  radar  is  functionally  represented  by  hierarchical  symbols  known  as  glyphs).  A 
radar  model  can  be  constructed  graphically  by  placing  and  connecting  glyphs  of  the  radar 
hardware,  environment  and  signal  processing  similar  in  principle  to  an  operational  radar 
system.  Not  limited  to  just  constructing  simple  radar  models,  the  functionality  available 
in  RLSTAP  allows  the  user  to  model  complex  radar  systems  in  defined  geographical  lo¬ 
cations  (site-specific  clutter) .  RLSTAP  is  a  time  based  simulation  in  which  the  positions 
of  the  platform  and  targets  change  from  pulse  to  pulse,  providing  excellent  simulation  of 
range  (slant  range)  and  Doppler  walk.  In  view  of  the  high  fidelity  simulation,  RLSTAP 
requires  significant  effort  in  understanding  radar  system  details  and  computation  resources 
to  develop  a  workable  and  realistic  radar  model. 

This  report  provides  the  RLSTAP  user  with  an  introduction  to  the  software  develop¬ 
ment  tools  currently  available  within  the  KHOROS  environment  to  facilitate  the  integra¬ 
tion  of  non-KHOROS  programs  into  the  KHOROS  environment.  The  KHOROS  environ¬ 
ment  has  a  complete  set  of  software  development  tools  for  modifying  existing  stand-alone 
program  written  in  either  “C”  or  “C++”  for  KHOROS  compliance.  Associated  manuals 
and  help  documents  are  automatically  generated  when  KHOROS  programs  are  created. 
These  documents  can  be  edited  manually  to  provide  additional  descriptions  and  examples 
to  the  user. 

We  have  modified  an  existing  stand-alone  modified  composite  (MC)  sea  clutter  pro¬ 
gram  written  in  “C”  for  KHOROS  compliance  and  integrated  the  model  with  existing 
RLSTAP  programs  with  minimum  code  modification.  This  exercise  has  allowed  us  to 
examine  the  ability  to  add  functionality  and  to  further  enhance  RLSTAP  modelling  ca¬ 
pabilities.  The  MC  sea  clutter  model  when  incorporated  into  the  KHOROS  environment 
provides  a  more  accurate  modelling  of  the  microwave  returns  from  the  sea  surface  for  radar 
operating  in  the  L-band  wavelengths. 
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Glossary 

0  grazing  angle  (to  horizontal  axis,  degrees) 

&i  incidence  angle  (to  vertical  axis,  degrees)  =  (90°  —  0) 

A  RF  wavelength  of  transmitted  signal  (m) 

<7°  mean  clutter  radar  cross-section  per  unit  area  of  illuminated  surface  ( dBm2m~ 2) 

/  frequency 

Cantata  A  visual  programming  environment  workspace  within  the  KHOROS  system. 
Composer  Editing  and  management  of  KHOROS  software  objects. 

Craftsman  Create  and  manage  KHOROS  toolbox  and  software  objects. 

CLUI  Command  Line  User  Interface 
DEM  Digital  Elevation  Model 

glyph  graphical  representation  of  program  modules  in  RLSTAP 

GUI  Graphical  User  Interface 

GUISE  Graphical  User  Interface  Specification  Editor 

LULC  Land  Use  Land  Cover 

LUT  Look-Up  Table 

MC  Modified  Composite  Sea  Clutter  Model 

KHOROS  A  software  integration  and  development  environment 

RLSTAP  Research  Laboratory  Space  Time  Adaptive  Processing 

TSC  Technology  Service  Corporation 

TTCP  The  Technical  Cooperation  Panel 
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1  Introduction 

This  work  was  performed  under  the  task  DST99/045- “RADAR  SYSTEMS  MOD¬ 
ELLING/PREDICTION  /VALIDATION”  for  TTCP  SEN-TP3  Task  4  Collaborative  Radar 
Simulation  Validation. 

Research  Laboratory  Space  Time  Adaptive  Processing  (RLSTAP)  is  a  compilation  of 
a  multiple-channel  radar  environment,  signal  and  hardware  programs  integrated  within 
the  KHOROS  environment.  These  programs  are  graphically  represented  as  glyphs  in  the 
KHOROS  environment.  A  radar  model  can  be  constructed  graphically  by  placing  and 
connecting  glyphs  of  the  radar  hardware,  environment  and  signal  processing  similar  in 
principle  to  an  operational  radar  system.  Not  limited  to  just  constructing  simple  radar 
models,  the  functionality  currently  available  in  RLSTAP  allows  the  user  to  model  complex 
radar  systems  in  defined  geographical  locations  (site-specific  clutter).  RLSTAP  is  a  time 
based  simulation  in  which  the  positions  of  the  platform  and  targets  change  from  pulse  to 
pulse,  providing  excellent  simulation  of  range  (slant  range)  and  Doppler  walk. 

This  report  discusses  an  approach  used  to  create  an  L-band  sea  clutter  model  in  the 
KHOROS  environment  and  integration  of  the  model  with  existing  RLSTAP  programs 
or  kroutines.  The  majority  of  programs  in  the  KHOROS  environment  are  kroutines. 
Currently,  RLSTAP  supports  two  programs  for  determining  clutter  scattering  from  the 
sea  surface:  (1)  the  “RSS  SigmaO  Curves”  program  generates  the  backscattering  coefficient, 
(j°  Look-Up  Table  (LUT)  using  the  Technology  Service  Corporation  (TSC)  RSS  model; 
and  (2)  the  “SigmaO  Curve”  program  generates  the  a°  LUT  using  either  the  ERAM  model 
or  the  Five  Parameter  model  or  a  user-specified  ASCII  Look-Up  Table. 

The  above  two  clutter  scattering  programs  are  representative  of  the  clutter  models 
commonly  used  in  the  radar  community.  They  are  however,  not  an  accurate  representation 
of  the  L-band  sea  clutter  characteristic.  We  have  developed  a  second-order  modified 
composite  sea  clutter  model  (MC),  which  is  based  upon  the  concept  of  a  two-scale  sea 
surface  model  to  predict  L-band  radar  sea  backscattering  coefficient,  a0,  for  varying  sea 
conditions,  grazing  angles  and  radar  look  directions.  We  have  compared  the  MC  sea 
clutter  model  and  commonly  used  semi-empirical  sea  clutter  models  with  data  collected  by 
the  Naval  Research  Laboratory  four-frequency  radar  system  under  varying  sea  conditions 
and  Jet  Propulsion  Laboratory  AirSAR  polarimetric  synthetic  aperture  radar  images  of 
the  North-West  of  Australia.  The  results  of  the  comparison  have  shown  that  the  MC  sea 
clutter  model  better  predicts  the  backscattering  coefficient,  a°,  of  the  sea  surface  returns 
for  a  wide  range  of  grazing  angles,  radar  look  directions  and  sea  conditions. 

The  first  two  sections  of  the  report  briefly  describe  the  RLSTAP  sea  clutter  programs 
and  the  formulation  of  the  second  order  L-band  MC  sea  clutter  model.  Section  4  of  the 
report  provides  the  reader  with  a  tutorial  like  approach  to  incorporating  new  glyphs  into 
the  KHOROS  environment.  Readers  not  interested  in  creating  new  glyphs  may  skip  the 
entire  section  and  proceed  to  the  modelling  results  described  in  Section  5.  The  approach 
described  in  this  report  is  by  no  means  the  only  method  available  for  integrating  new 
functionality  with  existing  RLSTAP  programs,  but  rather  a  practical  method  for  quickly 
converting  existing  stand-alone  “C”  programs  into  KHOROS  kroutines  or  glyphs. 
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2  RLSTAP  Clutter  Scattering  Models 


For  a  radar  system  operating  close  to  the  coastal  region,  or  inland  areas,  its  detection 
and  tracking  capabilities  will  be  severely  affected  by  compounded  microwave  return  signals 
from  the  coastal  water,  foreshore  vegetation,  land  cover  vegetation,  land  cover  usage  and 
man-made  structures.  The  clutter  returns  from  the  coastal  and  inland  areas  will  depend  on 
the  terrain  characteristics  e.g.,  vegetation  type,  presence  of  man-made  structures,  ground 
surface  roughness  and  many  more.  One  of  the  attractions  in  using  RLSTAP  for  radar 
modelling  is  the  availability  of  a  site-specific  scattering  module.  This  capability  not  only 
allows  for  more  accurate  modelling  of  the  microwave  returns  from  the  land  and  sea  surface, 
but  also  accounts  for  terrain  masking  or  shadowing. 

The  site-specific  clutter  modelling  module  allows  the  user  to  specify  the  area  of  interest 
by  defining  the  geo-location,  e.g.,  South-West  (reference  origin)  and  the  North-East  lati¬ 
tude  (-90,90)  and  longitude  (-180,180)  in  degrees.  Figure  2.1  shows  the  flow  process  of  the 
site-specific  clutter  modelling.  The  program  locates  and  extracts:  (1)  Digital  Elevation 
Model  (DEM)  data;  and  (2)  Land  Use  Land  Cover  (LULC)  information  for  the  specific 
region  of  interest,  and  processes  them. 


Site-Specific  $ 
Scatter  Field  Profile  j 

'  f 


Figure  2.1:  RLSTAP  Site  Specific  Clutter  Program  Flow 


There  are  however  a  few  constraints  associated  with  the  site-specific  module:  (1)  the 
resolution  of  the  DEM  and  LULC  data  is  at  present  fixed  at  3  arc-seconds  or  approximately 
90  m;  (2)  it  only  accepts  US  Geophysical  Service  (USGS)  DEM  and  LULC  data  formats; 
and  (3)  the  sea  and  land  clutter  scattering  models  are  limited.  The  accuracy  of  modelling 
the  microwave  returns  from  the  sea  surface  and  land  cover  is  dependent  on  the  availability 
of  a  suitable  scattering  model  in  RLSTAP. 

The  current  RLSTAP  version  3.13,  has  two  glyphs  or  programs  for  defining  clut¬ 
ter  scattering  from  the  sea  surface.  Both  glyphs  generate  a  LUT  of  the  backscattering 
coefficient,  cr°  versus  grazing  angle,  ip. 
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2.1  RSS  SigmaO  Curve 


The  Graphical  User  Interface  (GUI)  for  the  first  clutter  scattering  program  labelled 
“RSS  SIGMAO  Curve”  is  shown  in  Figure  2.2.  The  “RSS  SIGMAO  Curve”  glyph  generates 
a  LUT  of  the  backscattering  coefficient,  <r°  versus  grazing  angle,  ip,  for  a  specified  material 
code  (sea  or  land  cover  type)  using  the  TSC  RSS  model  for  grazing  angles  spanning  0  to  90 
degrees.  In  the  “RSS  SigmaO  Curve”  glyph,  the  user  has  the  option  to  generate  microwave 
returns  from  either  the  land  cover  or  from  the  sea  surface. 


For  land  clutter,  five  different  land  cover  types  are  available:  (1)  Desert /Roads;  (2) 
Rangeland;  (3)  Agriculture;  (4)  Forest;  and  (5)  Urban.  This  scattering  model  includes  the 
effect  of  seasonal  variation  on  the  strength  of  the  clutter  returns. 


For  sea  clutter,  six  different  sea  states  are  available,  ranging  from  sea  state  0  to  sea 
state  5. 


Figure  2.2:  GUI  for  RSS  SigmaO  Curve 
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2.2  SigmaO  Curve 

The  interface  to  the  second  clutter  scattering  program  labelled  “SIGMAO  Curve”  is 
shown  in  Figure  2.3.  The  “SIGMAO  Curve”  glyph  generates  a  LUT  of  a°  as  a  function  of 
grazing  angle,  x/>  (spanning  0  to  90  degrees),  using  one  of  the  following  clutter  models:  (1) 
ERAM,  (2)  Five  parameter  model;  or  (3)  user  specified  ASCII  LUT. 

In  the  ERAM  model,  the  user  has  the  option  to  generate  either  backscattering  from  the 
land  cover  or  from  the  sea  surface  for  a  radar  operating  in  either  X,  Ka  or  Ku  band.  The 
ERAM  model  was  provided  by  Raytheon  Systems  Company  without  any  documentation. 
It  is  not  suitable  for  modelling  radars  operating  at  frequencies  less  than  X-band. 

The  Five  parameter  model  is  independent  of  frequency,  with  the  backscattering  values 
determined  by 


<7°  =  cq  +  ci  sin^  +  C2exp{~(c3(90-V’)]C4} 


(1) 


where  co  =  Plateau,  c\—  Constant  Gamma  Multiplier;  C2  =  Exponential  Multiplier;  C3 
Incidence  Angle  Multiplier;  and  c\  =  Exponent. 


■  -  Generates  slgtnao  as  a  function  of  grazing  angle  based  on  one  of  the  following:  eram.  s  parameter  model,  or  use 


Ganaratat  slgmaD  at  a  Junction  of  grazjng  anfllt  basad  ononaof  tha  fotovrtng:  ERAM,  5  ptrarwtar  modal,  or  utar  tptcffltd  took- up  UMa  (ICON 


curva)  I 


Figure  2.3:  GUI  for  SigmaO  Curve  Model 

The  third  option  is  for  the  user  to  specify  an  ASCII  file  containing  a  list  of  grazing 
angles  followed  by  the  corresponding  a0  values. 
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3  L-band  Sea  Clutter  Model 

The  principles  of  microwave  scattering  by  slightly  rough  surfaces  have  been  applied 
to  backscattering  from  the  sea  surface  for  L-band  VV-polarisation  under  standard  at¬ 
mospheric  conditions  [2] .  Based  on  the,  two-scale  sea  surface  model  proposed  by  Guinard 
and  Daley  [2]  and  the  surface  scattering  model  proposed  by  Chen,  Fung  and  Weissman  [5], 
the  second-order  modified  composite  (MC)  sea  clutter  model  takes  into  account  polarisa¬ 
tion,  wind  speed,  azimuth  and  grazing  angles  variations.  In  this  model,  the  backscattering 
coefficient,  <7°,  can  be  represented  as  the  sum  of  separate  contributions  from  the  specular 
return  and  resonant  Bragg-scattering. 

Three  commonly  used  sea  clutter  models:  (1)  Georgia  Institute  of  Technology  (GIT): 
(2)  the  Hybrid  (HYB);  and  (3)  the  Technology  Service  Corporation  (TSC)  sea  clutter 
models  were  also  considered  for  comparison  [6].  These  sea  clutter  models  are  compared 
with  four  sets  of  measured  L-band  VV-polarised  data. 

Of  the  four  data  sets  used,  three  sets  of  calibrated  airborne  L-band  data  were  collected 
by  the  Naval  Research  Laboratory  (NRL)  in  the  late  60s  and  early  70s  from  several  sea 
locations  and  in  varying  sea  conditions.  The  fourth  set  of  data  was  obtained  from  mea¬ 
sured  L-band  VV-polarised  polarimetric  synthetic  aperture  radar  (POLSAR)  images  of 
the  Darwin  region.  These  images  were  acquired  during  the  PACRIM  mission  by  the  Jet 
Propulsion  Laboratory  on  the  23rd  November  1996. 


3.1  Data  Source  for  L  Band 

Four  sets  of  measured  L-band  VV-polarised  data  have  been  used  to  test  the  accuracy 
of  the  GIT,  HYB,  TSC  and  the  proposed  MC  sea  clutter  models.  Brief  descriptions  of 
the  NRL  4FR  data  sets,  and  the  DARWIN  POLSAR  data  are  summarised  below. 

3.1.1  Naval  Research  Laboratory  Airborne  Trials  [1,  2,  3,  4] 

The  Naval  Research  Laboratory  (NRL)  measurements  from  the  60s  and  early  70s 
used  for  validating  the  model  were  obtained  using  a  four-frequency  (4FR)  airborne  radar 
system.  The  carrier  frequencies  were:  (1)  UHF  (428  MHz);  (2)  L-band  (1.228  GHz);  (3) 
C-band  (4.455  GHz),  and  (4)  X-band  (8.910  GHz) [7,  3,  4],  For  all  analyses,  the  median 
sea  backscattered  values  were  converted  to  mean  values  assuming  a  Rayleigh  amplitude 
distribution  [7]. 

1.  NRL-4FR  Data  [1] 

Sea  backscattered  measurements  were  recorded  using  the  NRL  4FR  radar  system, 
off  New  Jersey  in  1964,  off  Puerto  Rico  in  1965  and  in  the  North  Atlantic  in  1969. 
Surface  truth  information  for  the  Puerto  Rico  deployment  was  measured  by  a  team 
of  researchers  from  the  Applied  Physics  Laboratory,  Johns  Hopkins  University,  in 
an  instrumented  vessel  [2].  For  the  North  Atlantic  deployment,  the  surface  truth 
information  was  provided  by  oceanographic  surface  vessels  operated  by  the  French 
and  British  governments  [2],  We  will  hereafter  refer  to  data  measured  as  the  NRL- 
4FR  data. 
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2.  Joint  Ocean  Surface  Study,  1970  (JOSS  I)  [4] 

The  NRL  4FR  system  was  used  to  measure  sea  backscattered  signals  off  the  vicinity 
of  Argus  Island,  Bermuda.  Wind  and  wave  height  information  were  recorded  by 
observers  stationed  on  Argus  Island. 

3.  Joint  Ocean  Surface  Study,  1971  (JOSS  II)  [8] 

The  NRL  4FR  system  was  used  to  measure  sea  backscattered  signals  off  the  East 
Coast  of  Northern  America.  Wind  velocity,  direction,  and  average  wave  height 
information  were  obtained  hourly  from  an  instrumented  buoy  and  an  ocean  station 
vessel.  [8j. 


3.1.2  AirSAR  Polarimetric  Synthetic  Aperture  (POLSAR)  Data 

POLSAR  is  a  NASA-sponsored  polarimetric  synthetic  aperture  radar  system  developed 
at  the  Jet  Propulsion  Laboratory  in  California.  The  POLSAR  system  provides  high  quality 
polarimetric  SAR  data  in  three  frequencies;  P  (0.44  GHz),  L  (1.25  GHz)  and  C  (5.39  GHz) 
bands  in  four  transmit/receive  polarisations  (VV,  HH,  VH,  HV). 

The  study  area  of  interest  was  in  the  vicinity  of  Darwin  (130°  45  ’  E,  12  °  30’  S)  in 
the  Northern  Territory  of  Australia.  Two  flight  passes  were  selected  and  the  two  L-band 
polarimetric  SAR  images  of  vertical  (VV)  polarisation  collected  with  a  bandwidth  of  40 
MHz  are  shown  in  Figures  3.1(a)  and  (b).  Specifications  of  the  POLSAR  images  are 
tabulated  in  Table  3.1.  The  number  of  looks  processed  in  azimuth  is  18,  and  1  in  range. 
Recorded  measurements  of  the  sea  conditions  are  tabulated  in  Table  3.2. 


Table  3.1:  POLSAR  Specifications 


Description 

CM5512 

CM5148 

Near  Slant  Range  (m) 

9013.26 

9013.26 

Far  Slant  Range  (m) 

17537.36 

19244.51 

Near  Look  Angle  (0*) 

22.5° 

22.6° 

Far  Look  Angle  ( ) 

61.7° 

61.7° 

Range  Resolution  (m) 

3.3310 

3.3310 

Azimuth  Resolution  (m) 

9.2592 

9.2592 

Table  3.2:  Hydro-meteor  Parameters 


Description 

Darwin  Harbour 

Fannie  Bay 

Air  Temperature  (Dry) 

35.0°C 

32.5°C 

Air  Temperature  (Wet) 

29.0°C 

27.0°C 

Relative  Humidity 

68% 

64% 

Wind  Speed 

14  km/h 

4  kt 

Wind  Direction 

O 

O 

IV 

y*H 

SE 

Wave  Height 

1  cm 

Calm 
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(b)  CM5148 

Figure  3. 1:  POLSAR  Darwin  L-band  VV  data 
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3.2  Modified  Composite  (MC)  Sea  Clutter  Model 

The  MC  sea  clutter  model  is  based  upon  the  concept  of  slightly  rough  sea  surface 
scattering  in  which  the  dominant  mechanism  is  Bragg  resonant  scattering,  whose  am¬ 
plitude  is  related  to  the  second-order  directional  sea  spectrum  proposed  by  Chen,  Fung 
and  Weissman  [5].  The  model  also  accounts  for  specular  return  near-normal  incidence. 
The  backscattering  coefficient,  cr°,  from  the  sea  surface  can  be  represented  as  the  sum  of 
separate  contributions  from  specular  return  and  Bragg-scattering. 

a°  =  (2) 

Specular  Returns  Bragg  Scattering 


The  MC  sea  clutter  model  is  constructed  based  on  these  two  forms  of  scattering  mecha¬ 
nism.  The  primary  scattering  mechanism  in  the  low-grazing  angle  and  plateau  regions  is 
assumed  to  be  Bragg  scattering  whereby  the  backscattering  coefficient,  a°B,  is  in  propor¬ 
tion  to  the  spectral  density  of  the  short  gravity-capillary  waves  for  a  fully  developed  sea 
[9].  At  near  normal  incidence,  the  effect  of  specular  return,  a °s,  due  to  facets  normal  to 
the  radar  signal  is  included. 


3.2.1  Specular  Backscattering  for  Fully  Developed  Sea 

Near-normal  incidence  backscatter  from  the  sea  surface  is  better  described  by  the 
specular  point  or  tilted  facet  model.  For  an  isotropic  rough  surface  of  Gaussian  statistics, 
the  backscatter  return  at  near  normal  incidence  is  given  by  [10] 

cr%  =  sec4  &i  exp(—  tan2  6i/s2)  (3) 

sz 

where  R(0)  is  the  Fresnel  reflection  coefficient  of  the  sea  surface  at  normal  incidence,  s2 
is  the  mean-square  sea  slope  and  8i  is  the  angle  of  incidence. 

(a)  Mean-Square  Sea  Slope,  s2 

2  f  5.12  x  10 ~3U,  U  <  5m/s 
S  ”  1 1.56  x  10 -3U,  U  >  5m/s 

where  U  (m/s)  is  the  wind  speed  at  an  altitude  12.5  m  above  the  sea  surface. 

(b)  Reflection  Coefficient,  |R(0)| 


\m\  = 


0.52 


er  —  1 


where  er  is  the  relative  dielectric  constant  of  sea  water. 


(5) 
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3.2.2  Bragg  Scattering 

1.  Second  Order  Bragg  Scattering  Model, 

Horizontal  Polarisation,  aB2(HH)  =  ko  sin4(90  -  9i)aHH  W{K,<f>)  (6) 


Vertical  Polarisation,  a 


B2{W) 


Directional  Sea  Spectrum 

k20  sin4 (90  -  9i)  aVv  W(K,<j>) 

Directional  Sea  Spectrum 


where  k0  =  ^  is  the  radar  wavenumber,  9{  is  the  angle  of  incidence,  K  is  the  sea 
wavenumber,  and  <j>  is  the  wind  direction  with  respect  to  boresight. 

2.  The  first-order  scattering  coefficient,  app  is  [2] 

(er  —  1) 


Horizontal  Polarisation,  oihh  = 
Vertical  Polarisation,  ayv  — 


[cos^i  +  \/(er  -  sin2#*)]2 
(er  —  l)[er(l  +  sin2  0i)  —  sin2  9{ 
[er  cos  9{  +  y/(er  -  sin2#*)]2 


(7) 


where  er  is  the  relative  dielectric  constant  of  the  sea  surface  with  value  72  —  j59 
(L-band,  1.275  GHz)  [9]. 

3.  Two-dimensional  (2-D)  Sea  Ripple  Spectrum 

Adopting  the  2-D  sea  ripple  spectrum  developed  by  Pierson  and  modified  by  Chen 
and  Fung  [5],  the  2-D  directional  (taking  into  account  the  wind  direction)  sea  spec¬ 
trum  is 


W{K,<f>)  =  —  (l  +  rcos20)W(AT) 
2n 


(8) 


(a)  The  constant,  r,  is  given  by 


(IT*') 

where  v  is  the  ratio  of  the  slope  variances  in  the  crosswind  and  upwind  directions 

[5], 


v  = 


0.003+ 1. 92 x10~3E7  jt  ^  c  / 
3.16x10-3  U  ’  U  ^  om/  o 


3.16xl0-3£7 
0.003+0.84  xl0~317 


0.003+0.84  X  IQ- J  (7  tj  >  r  / 
0.005+0. 78  xlO-3C/’  U  —  jm/& 


(b)  Wind  Direction,  0 


{0°,  upwind 

90°or  —  90°,  crosswind 
180°,  downwind 

(c)  1-D  sea  ripple  spectra  for  a  single  wave  train,  W{K) 


W(K)  =  A(2tt)p 


-i  (1  +  3&9LYl 


rQ-922  j R  ^  +  ig.) 


(i+p) 


(9) 
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where  A  =  1.75  if  U  <  5  m/s  or  A  =  0.875  if  U  >  5  m/s,  K ^  =  gp/r  = 
(3.63)2cm~2,  g  =  acceleration  of  gravity  (9.81  m/s2),  p  =  sea  water  density,  and 
r  =  sea  surface  tension.  Note  that  according  to  Valenzuela  [11],  K  should  have 
a  —3.922  power  law  rather  than  the  usual  -4  power  law  for  L-band  resonant 
Bragg  waves. 

The  relationship  between  the  sea  surface  wind  speed  at  an  altitude  z  and  friction 
velocity  at  neutral  stability  of  the  atmosphere,  U*  is  [5] 


‘'-S'-i 


where 

Z0(cm)  = 

and  p  =  5  —  log10  U* . 


0.684 

U* 


+  4.28  x  10 ~5C/*2  -  0.0443 


3.3  Comparisons  with  Radar  Measurements 

Plots  of  the  results  in  Figures  3.2  (a)  and  (b)  show  good  agreement  between  the 
measured  NRL-4FR  data  and  the  MC  predicted  <r°  values  for  the  two  wind  speeds.  Both 
the  HYB  and  TSC  models  underestimate  the  backscattered  values,  a0. 

We  further  tested  the  validity  of  the  MC  sea  clutter  model  using  the  JOSS  I,  JOSS  II, 
and  the  POLSAR  DARWIN  measured  data.  Figure  3.3(a)  again  shows  the  poor  agreement 
of  the  TSC  and  HYB  models  with  the  JOSS  I  measured  data.  Where  these  two  models 
fail,  the  MC  sea  clutter  model  shows  good  agreement  with  the  measured  data.  Figure 
3.3(b)  shows  that  the  MC  sea  clutter  model  accurately  models  the  JOSS  II  measured 
data.  Note  that  the  GIT  model  has  good  agreement  with  the  NRL-4FR  and  JOSS  I  data 
for  limited  grazing  angles,  €E  [5°,  10°].  However,  further  validation  is  required  due  to  the 
lack  of  low  grazing  angle  data  in  the  available  data  sets. 

For  the  POLSAR  data  acquired  over  the  DARWIN  region  under  low  wind  speed  condi¬ 
tion,  the  MC  sea  clutter  model  accurately  predicted  the  backscattering  coefficient,  a°,  for 
20°  <  4>  <  65°.  Figure  3.3(c)  shows  the  measured  data,  and  the  predicted  a°  values  by  the 
MC  and  semi-empirical  sea  clutter  models.  The  TSC  and  HYB  models  underestimate 
the  backscatter  values  when  compared  to  the  measured  data. 

We  have  tested  four  sea  clutter  models  against  measured  data.  These  results  suggested 
that  among  these  four  (GIT,  HYB,  TSC,  MC)  sea  clutter  models,  the  MC  sea  clutter 
model  is  the  preferred  generic  sea  clutter  model  for  predicting  the  L-band,  VV-polarised 
sea  clutter  for  a  wide  range  of  grazing  angles. 
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4  L-band  MC  Sea  Clutter  Glyph  Creation 

This  section  illustrates  an  application  of  the  KHOROS  software  development  tools  to 
create  a  new  sea  clutter  toolbox  to  contain  the  MC  sea  clutter  model  and  a  method  to 
convert  an  existing  MC  sea  clutter  program,  labelled  “SeaClutter.c”  into  a  KHOROS 
program  or  kroutine.  In  the  KHOROS  system,  a  toolbox  is  a  collection  of  programs  and 
libraries  grouped  together  with  a  common  functionality.  The  majority  of  these  programs 
are  in  the  standard  KHOROS  kroutines  format. 


Step  1.  Toolbox  Creation 

The  first  step  in  inserting  a  “stand-alone  C”  code  into  the  KHOROS  environment 
is  to  create  either  a  new  toolbox  or  use  an  existing  toolbox.  This  is  done  using 
Craftsman  [12,  13].  The  Craftsman  application  is  used  to  create  and  manage 
KHOROS  toolbox  and  software  objects.  Start  up  Craftsman  from  the  command 
line.  When  Craftsman  is  displayed  (refer  to  Figure  4.1),  you  will  notice  two  lists: 
(1)  the  list  on  the  left  corresponds  to  the  toolboxes  currently  available;  and  (2)  the 
right  list  corresponds  to  a  selected  toolbox  software  objects. 


Crsfararv'ioftv/c  re  Development  EnvIrcFT1! 


In*.*.  v|  VMfeltan  (ggg  m 

Toolboxes  <JJ>  f  (no  toolbox  selected)  f 


e.ebout 

a_  ant  on  a_  server,  1  Ibran 
a,con  trailer*  ■. 
4_d»te_Mnlp 
a.dwlgn.ejr'thesu 


a.d  1  s  1  ta  l  .da  taj#*  oce&s  1 1 
e.dlfl ital.al 90a l _proc*« 
a.env  1  ron«nt.saner  aton 


Toolbox  Operations  Vf 

Create  loci bo*  J 

..Adel  Toolbox  Reference 

lCrA4t4:  idclt»i-  niniiil 
K&iM,  Toole*-  Rif’d/urstp.  . 


Figure  4- 1-  The  Craftsman  software  development  environment 


(a)  To  create  a  new  sea  clutter  toolbox,  select  Create  Toolbox  from  the  Toolbox 
Operations  pulldown  list. 

(b)  This  action  brings  up  a  sub-menu  shown  in  Figure  4.2  and  allows  you  to  enter 
all  the  software  object  information. 

1.  Toolbox  Name 

Type  in  the  toolbox  name  as  you  would  like  it  to  appear  in  Craftsman. 
In  this  example,  the  sea  clutter  toolbox  will  appear  on  the  Craftsman 
toolboxes  list  as  MC. 
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Figure  4-2:  The  Craftsman  Create  Toolbox  Sub-menu 


2.  Toolbox  Path 

Directory  preference  for  toolbox  location. 

3.  Toolbox  Title 

Short  description  of  the  toolbox  title. 

4.  Point  of  Contact  Information 

The  point  of  contact  information  contains  the  author’s  name  and  email  ad¬ 
dress,  primarily  used  for  bug  control  reporting  for  public  release  KHOROS 
programs. 

(c)  Click  on  the  Create  Toolbox  button,  and  the  new  sea  clutter  toolbox  labelled 
MC  is  now  created.  It  will  appear  in  the  toolbox  list  as  shown  in  Figure  4.3. 


i-lcraftsman:  Software  Development  Envlrcl 


Options  vj 

Toolboxes  <34>  f 


,turj 


f  Help  l|  Quit 


MC  (all:  0)  f 


’envision" 

geonetry 

lnage 

eetrlx 

■ 

meat,,; 

migration 
rotro v  .V' 

AMolwlAhA 


- — — - 1 

L 

Toolbox  Operation*  Vj 

Object  Operations 

A 

Figure  4-3:  The  Newly  created  Sea  Clutter  Toolbox  labelled  MC 
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Step  2.  Progam  Object  Creation 

The  next  step  involves  creating  the  sea  clutter  KHOROS  program  object.  The 
program  object  comprises  the  MC  sea  clutter  program  along  with  a  KHOROS 
imposed  directory  structure  and  KHOROS  generated  code. 

(a)  Select  Create  Object  from  the  Object  Operations  pulldown  list.  Selecting 
this  will  bring  up  the  object  creation  submenu  as  shown  in  Figure  4.4(b). 


(a) 


(b) 


Figure  4-4'-  (a)  Craftsman  ;  and  (b)  the  Object  Creation  Submenu 


(b)  Fill  in  all  the  program  objects  attributes. 

1.  Object  Name 

Title  of  the  object  and  the  source  code  of  the  MC  sea  clutter  program  will 
be  placed  in  a  directory  of  this  name. 

2.  Binary  Name 

Name  of  the  executable  and  library  archive  (if  not  specified,  default  to 
Object  Name) . 
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3.  Icon  Name 

Title  of  glyph  which  appears  on  the  Cantata  pulldown  list  as  a  subcategory 
(if  not  specified,  default  to  Object  Name). 

4.  Category  and  Subcategory 

Use  to  classify  the  object  to  a  broad  domain  application.  In  this  example, 
the  MC  sea  clutter  model  is  classified  under  a  broader  category  of  “Sea 
Clutter” .  Within  the  Sea  Clutter  category,  the  MC  sea  clutter  model  falls 
into  the  subcategory  of  “L-band  Sea  Clutter”. 

5.  Generated  Language 

Choices  are  “C”  or  “C++”.  The  existing  sea  clutter  program  was  written 
in  “C”,  therefore  the  “C”  option  was  selected. 

6.  Install  in  Cantata 

This  toggle  is  used  to  specify  whether  you  want  the  object  to  be  visible  to 

Cantata  . 

7.  Create  Library  Routine 

This  toggle  is  used  to  specify  whether  you  want  a  library  object  to  be 
associated  with  your  object. 

Upon  completion,  click  on  the  Create  KROUTINE  to  create  the  MC  sea 
clutter  program  object  (directory  structure,  makefiles,  source  code  and  others 
associated  with  the  program  object).  At  this  stage,  the  main  KHOROS  source 
file  generated  for  the  MC  sea  clutter  kroutine  is  labelled  “MC.c”  and  the  object 
code,  “MC.o”.  When  the  new  MC  program  object  has  been  created,  it  will 
appear  highlighted  in  the  Craftsman  object  list  as  shown  in  Figure  4.5. 


Craftsman:  Software  Development  Envirc 


Help  ■  Quit 


Options  VI 


Toolboxes  (34)  f 


HC  (all:  1)  f 


image 

matrix  -  . 

mean 

migration 

retro 

sampledata 

• 

support 

Toolbox  Operations 


1 


Object  Operations 


3 


Figure  4-5:  Program  Object  Creation 
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Step  3.  Create  and  Edit  Graphical  User  Interface 

When  Craftsman  creates  the  MC  program  object,  it  copies  a  generic  or  template 
User  Interface  Specification  (UIS)  “MC.pane”  file  to  the  new  program  object  di¬ 
rectory.  The  “MC.pane”  file  specify  the  input  and  output  arguments  of  the  MC 
program.  Since  the  user  interface  is  not  what  we  wanted,  we  need  to  edit  the 
“MC.pane”  using  either  a  text  editor  or  graphically  using  a  KHOROS  tool,  the 
Composer  application  to  provide  the  require  input  and  output  arguments. 

(a)  Select  the  Edit  Object  (Composer)  from  the  Object  Operations  pulldown 
menu  shown  in  Figure  4.6  (a). 


(a)  (b) 

Figure  4.6:  (a)  Edit  Object  (Composer)  and  (b)  Composer  Software  Object  Editor 


(b)  This  action  invokes  the  Composer  object  editor  shown  in  Figure  4.6(b).  Com¬ 
poser  editor  allows  editing  and  management  of  existing  software  objects.  In 
addition,  Composer  allows  GUI  editing  via  GUISE  (Graphical  User  Inter¬ 
face  Specification  Editor),  sets  software  objects  attributes,  generate  code,  and 
compiles  the  program. 

(c)  We  wish  to  modify  the  UIS  file  of  the  MC  program  interactively,  so  select 
UIS  under  the  List  Files  (high- lighted  by  default),  then  select  Guise  from 
the  File  Operations  drop-down  menu  as  shown  in  Figure  4.6(b).  GUISE  is 
a  graphical  design  tool  used  to  create  and  modify  user  interfaces  of  software 
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Graphical  user  interface  Specification  Editor  CGUISEJ 
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Figure  4-7:  GUISE  GUI  Design  Tool 


objects  interactively.  This  action  invokes  two  windows:  (1)  GUISE  editing  tools 
window;  and  (2)  a  standard  GUI  of  the  MC  program  is  shown  in  Figure  4.7 
(bottom  figure)  for  further  modification. 


Generates  L-band  Sea  Clutter  Sigmao  vs  Grazing  Angles  EUT 


Generates  L-band  Sea  Clutter  SignaO  vs  Grazing  ftigle 


Figure  If.. 8:  MC  Sea  Clutter  GUI 


(d)  After  some  editing,  the  final  MC  sea  clutter  GUI  is  shown  in  Figure  4.8. 

(e)  Whenever  the  UIS  file  or  user  interface  of  a  program  is  modified,  we  have 
to  regenerate  the  outdated  KHOROS  source  code  and  documentation.  Select 
Commands  from  the  right  side  of  Composer  to  invoke  the  Commands 
window.  Then  click  the  Generate  Code  button  to  update  all  source  code 
and  documentation  based  on  the  new  user  interface  as  specified  in  “MC.pane”. 
While  Composer  is  regenerating  all  source  code  and  documentation  for  the 
program,  it  will  display  the  regenerating  status  as  shown  in  Figure  4.9. 


18 


DSTO-TR-1206 


Figure  4.9:  Regenerating  MC  Sea  Clutter  GUI  code  and  documentation 

Step  4.  Edit  Source  Code 

At  this  stage,  we  can  run  the  MC  sea  clutter  glyph  in  Cantata  minus  the  function¬ 
ality.  Cantata  is  a  visual  programming  environment  within  the  KHOROS  system 
where  programs  are  represented  as  icons  or  glyphs.  Next,  we  must  link  the  main 
KHOROS  file,  “MC.c”,  to  our  sea  clutter  program,  “SeaClutter.c”,  and  modify  “Sea- 
Clutter. c”  for  KHOROS  compliance. 

(a)  Insert  the  function  prototype  for  sea  clutter  program, 
void  SeaClutter(void) ; 

in  the  main  KHOROS  file  “MC.c”.  Refer  to  Appendix  A. I  for  an  example. 

(b)  Insert  call  to  the  sea  clutter  program  in  the  main  KHOROS  file  “MC.c” 

In  the  main  KHOROS  source  file  labelled  “MC.c”,  insert  the  call  to  the  sea 
clutter  program  as  follows, 

/*  main_library_call  */ 

SeaClutter(); 

/*  main_library_call_end  */ 

Refer  to  Appendix  A.l  for  an  example. 

(c)  Insert  the  KHOROS  header  file, 

#include  "MC.h" 

in  “SeaClutter.c”,  so  that  the  routine  will  understand  KHOROS  function  calls.. 
An  example  of  a  KHOROS  generated  header  file  labelled  “MC.h”  can  be  found 
in  Appendix  A. 2. 
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(d)  Change  main()  in  “SeaClutter.c”  to  a  sub-routine  (Refer  to  Appendix  A. 3  for 
an  example). 

(e)  Change  file  declaration  in  “SeaClutter.c”  to  KHOROS  file  declaration  (Refer 
to  Appendix  A. 3  for  an  example). 

(f)  Comment  or  remove  any  default  values.  These  values  are  now  accessible  through 
the  MC  GUI. 

(g)  Reassign  KHOROS  Common  Line  User  Interface  (CLUI)  input  variables  to 
program  variables.  Note  that  CLUI  variables  are  declared  global  variables  and 
do  not  have  to  be  passed  during  function  calls. 

(h)  Change  file  I/O  in  “SeaClutter.c”  to  KHOROS  file  I/O. 

(i)  Edit  the  “Imakefile”  and  “Makefile”  to  include  both  the  MC  sea  clutter  source 
file  “SeaClutter.c”  and  object  file,  “SeaClutter.o” 

(j)  After  completing  all  the  above  editing,  we  need  to  recompile  the  source  code  and 
update  all  documentation.  In  the  Composer  window  click  on  the  Commands 
button  and  select  the  Make  Install  option.  After  compilation  (without  errors), 
the  KHOROS  MC  sea  clutter  program  can  now  be  run  in  Cantata. 

Step  5.  Documentation  of  the  MC  Sea  Clutter  Program 

Adequate  documentation  of  newly  created  or  updated  programs  is  an  important  as¬ 
pect  of  a  software  development  process  or  life-cycle.  The  KHOROS  system  provides 
a  framework  for  supporting  and  writing  documentation. 

During  the  creation  of  the  MC  software  object,  the  KHOROS  system  automatically 
generates  3  types  of  documentation  that  conform  to  the  KHOROS  documentation 
standard:  (1)  an  on-line  manual  or  man  page  similar  to  the  UNIX  man  page,  de¬ 
scribing  the  functionality  and  usage  of  the  MC  program;  (2)  an  on-line  help  page 
accessible  from  within  Cantata;  and  (3)  a  printed  manual  providing  hardcopy  pre¬ 
sentation  of  the  documentation. 

The  above  mentioned  documentations  can  also  be  edited  to  provide  additional  de¬ 
scriptions  and  examples  to  the  user.  KHOROS  has  two  categories  of  man  pages:  (1) 
one  for  programs  (“MC.l”  page);  and  (2)  one  for  libraries  (“MC.3”  page). 

(a)  Editing  MC  Sea  Clutter  Program  Man  Page 

All  man  pages  are  accessible  from  Composer  by  selecting  the  DOC  files  op¬ 
tion.  To  edit  the  “MC.l”  man  page,  select  “MC.l”  from  the  list,  then  select 
Edit  from  the  File  Operations  pulldown  list.  This  action  invokes  an  editing 
window  containing  the  “MC.l”  man  page  with  KHOROS  tags  that  denote  text 
segment.  Edit  the  “MC.l”  man  page  to  add  the  desired  documentation,  then 
save  the  update  upon  completion. 

(b)  Updating  MC  Help  Pages 

To  propagate  the  changes  made  in  the  “MC.l”  man  page  to  the  help  page, 
select  the  Generate  Code  from  the  Composer  Commands  button.  You 
are  able  to  view  the  newly  formatted  help  page  by  selecting  “MC.hlp”  from 
the  file  list  then  select  the  Formatted  option.  Figure  4.10  shows  an  updated 
“MC.hlp”  page  with  addition  documentation  highlighted  in  a  box. 
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1  PROGRAM 

i"  y,  k  so  ^  -oy  '  s  is 

MC  -  Generates  L-band  Sea  Ciutter  SfgmaO  vs  Grazing  Angles  LUT 

_  >  *».  if  ^  f)  A  „**  ”  ^  1*2)  >''-;/t ' 5 1.-  r~  ■.  '  <  ')  '  /' 

2.  DESCRIPTION  ;  '  ■  t  . 

The  MC  Sea-xlutter  model  Is  based  upon  the  concept  of  the  two-scale 
sea  surface  model  to  predict  the  L-band  radar  sea  surface  backseat-  , 
termg  coefficient,  SlgmaO.  ,<*  ,  3  *  ,  ,  «  *  ”  v 

3,  PANE  ;-,RGLMl:N’S 


'Transmit  Centre  Frequency"  (required  float,  f  >  0  0) 
float  'i  ;T  >..V  ^  ^ 


"POLARISATION"  (required  toggle). 

io),.  \ 

or  2  (2) 

Polarisation  of  co- polarised  transmit  radar  signal  (VV  or  HH)  ■ 
"WIND  DIRECTION"  (required  toggle) 

mo).' 

2  (90), 
or  3  (1 80) 

Integer  toggle  selection 

"Wind  Speed"  (required  float,  WindSpeed.>=  C.Q) 

Measured  wind  speed.  above  sea  surface  in  m/s  , 

"Sea  Ciutter  Output"  (required  outfile) 

SlgmaO  versus  grazing  angles  :..J. 

4.  EXAMPLES 

5.  SEE  ALSO 

6  RESTRICTIONS^  |§  189^8 

p  ' '  1  "S  t  -f  v  t  ■”l4  *■'  a  '■'}  t  1  *  '  \  *  „<■  y  >  v 

MC  Glyph  Copyright  (C)  2001  j  SSR,  SSD,  DSTO  Salisbury,  Australia. 

f'N.  ’  5$  v  1^,  1  f  L  4?  ^  S  ”'■»  4  l'  W||  s  * 

7.  REFERENCES- 

,  ■  <  •  ...... 

J,  Pop.Lian  Choong,,  "Modelling  Airborne  Radar.L-band  Sea  and  Coastal 


copyright  (Cj;  1 993  -  1 997,:  KhprajRf^arch,  Inc.  (”KRI");  Altrigjits 


Figure  4-10:  An  Example  of  the  MC  Help  Page 
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We  have  in  this  section  describe  a  step-by-step  approach  to  modifying  and  incor¬ 
porating  an  existing  stand-alone  non-KHOROS  program  into  the  KHOROS  envi¬ 
ronment.  The  process  of  generating  man  and  help  pages  for  associated  KHOROS 
programs  have  also  been  briefly  discussed.  Figure  4.11  summarises  the  process  of 
incorporating  user  stand-alone  “C”  code  into  the  KHOROS  environment. 


In  Composer,  select  the  "DOC"  files 
option,  Then  select  Edit  from  File 
Operations  pulldown  menu 
to  update  documentatio 


Run  Craftsman  Software  Development  Environment 
Then  select  Create  Toolbox  from  the  Toolbox 
Step  1  M  Operations  pulldown  list. 


Edit  KHOROS  generated 
"C  code  and  own  routine. 


In  Craftsman,  select 
Create  Object  from  the 
Objects  Operations 
pulldown  list. 


In  Craftsman,  select  Edit  Object 
( Compose r)  from  the  Objects 
Operations  pulldown  list.  Select 
GUISE  from  the  File  Operations 
pulldown  list. 


Figure  4- IF  Flow  Diagram  of  the  Glyph  Creation  Process 
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5  Running  the  MC  Sea  Clutter  Program  in 

Cantata 

Now  that  we  have  compiled,  installed  and  edited  the  man  and  help  pages  for  the 
MC  sea  clutter  program,  we  can  execute  the  MC  program  from  within  Cantata.  Exe¬ 
cute  Cantata  from  the  UNIX  command  line  and  wait  for  it  to  initialise.  The  Cantata 
workspace  consists  of  a  large  canvas  on  which  the  program  line-up  can  be  constructed. 
The  following  steps  illustrate  a  method  for  creating  a  program  line-up  in  Cantata  . 

Step  1.  Selecting  Glyph  for  Program  Line-Up 

Select  Glyphs  from  the  selection  bar  at  the  top  of  the  workspace.  Then  select  Sea 
Clutter  — >  L-band  Sea  Clutter  — >  MC  (category  and  sub-category  attributes 
that  we  have  previously  inserted  during  the  sea  clutter  toolbox  creation  process). 


Figure  5.1:  MC  Glyph  and  GUI  in  Cantata 


Step  2.  Positioning  the  Selected  Glyph  onto  the  Cantata  Workspace 

Highlight  MC  to  select  and  drag  the  MC  glyph  to  the  Cantata  workspace.  Click 
once  with  the  left  mouse  button  to  place  the  MC  glyph  at  the  desired  position.  You 
can  reposition  the  glyph  by  selecting  it  with  the  left  mouse  button  and  moving  it 
to  the  desired  position.  Release  the  mouse  button  to  place  the  glyph  in  the  new 
position. 

Step  3.  Pane  Information 

Figure  5.1  shows  the  MC  glyph  and  if  we  click  on  the  triangle  located  at  the  top- left 
side  of  the  glyph  (pane  access  button),  we  can  view  the  MC  GUI  and  all  default 
values.  These  values  can  be  updated  before  running  the  program. 

Click  on  the  Help  button  to  access  on-line  information  about  the  glyph  functionality 
and  variables  descriptions.  Figure  4.10  shows  the  “MC.hlp”  page  accessible  via  the 
MC  glyph  Help  button. 
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A  simple  program  line-up  to  plot  the  output  of  the  MC  glyph  is  shown  in  Figure  5.2. 
Click  on  either  the  run  button  located  at  the  centre  of  the  MC  glyph  to  execute  the  line-up 


|I  ;':hv:  :  [  IWMWIijpn:  j  | 


WCBIN/HC  yf  l.20**$9  -polar  Inti  on  2  -UlndBic  1  -UlndSpnd  3  -OutputIWu  Avm/choorqp/%j*/  10*580.0 

•  (ntw)  ftJWOOaOCY  ttflwQ.curv*  SI  CHAO  Ctrvo 

S^LPHDOCWUlCTlIN/aiBnCl.ourv*  -dug  1  -girt  1  -o.tlg  /faM/ctaongpA«p/lo4SlD.l  -tm5  -Ijna  /tan/ctworgp/UpAi 

•  (nw)  MlSU«.iaTI»  plotter  P»okt4f 

M.VISUa.IZfiTIOKilN/plotUr  -dU*-o4ticJ*v«l  1  -l  /htWctxo^/t-p/ro45M.l 


Figure  5.2:  Simple  MC  line-up 


or  the  RUN  button  from  the  selection  bar  at  the  top  of  the  Cantata  workspace.  While 
executing,  the  simulation  parameters  are  displayed  in  the  bottom  Cantata  workspace  as 
shown  in  Figure  5.2.  An  output  plot  of  the  backscatter  coefficient,  a°  versus  grazing  angle, 
ip  is  shown  in  Figure  5.3. 


Figure  5.3:  MC  Output  plot 
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5.1  A  Comparison  between  the  MC  Sea  Clutter  Model  and 
“RSS  SigmaO  Curve”  Model 

We  briefly  compared  the  output  of  the  MC  sea  clutter  model  and  the  “RSS  SigmaO 
Curve”  Model  for  modelling  the  backscattering  characteristics  from  the  sea  surface.  We 
choose  to  compare  these  models  for  two  sea  conditions:  (1)  Low  wind  speed  [wind  velocity, 
U  =  1.3  m/s  (sea  state  1)];  and  (2)  moderate  wind  speed  with  [wind  velocity,  U  =  5  m/s 
(sea  state  2)]. 

The  plots  shown  in  Figures  5.4  and  5.5  describe  the  a°  variation  as  ip  varies  from  low 
to  high  values.  The  “RSS  SigmaO  Curve”  sea  clutter  model  has  a  number  of  drawbacks: 
(1)  the  model  does  not  account  for  specular  returns  for  ip  at  near-normal  incidence  angles; 
and  (2)  the  backscatter  returns  are  not  representative  of  measured  L-band  VV-polarised 
airborne  radar  data  as  shown  in  Figures  3.2  (a)  -  (b)  and  Figures  3.3  (a)  -  (c).  . 


6  Conclusions 

This  report  focused  on  two  tasks:  (1)  the  modelling  of  L-band  VV-polarised  sea  clutter 
data;  and  (2)  the  ability  to  add  functionality  into  the  KHOROS  environment  for  integra¬ 
tion  with  RLSTAP  programs. 

In  modelling  L-band  sea  clutter  data,  the  results  of  the  comparisons  have  shown  that 
the  two-scale  model  modified  for  anisotropic  sea  spectrum,  the  MC  sea  clutter  model, 
accurately  predicted  the  NRL-4FR,  JOSS  I,  and  JOSS  II  data  for  a  wide  range  of  grazing 
angles  and  sea  conditions.  In  addition,  the  MC  sea  clutter  model  accurately  predicted 
the  measured  backscattering  coefficient,  a°,  for  the  sea  surface  around  the  Port  Darwin 
and  Fannie  Bay  areas  for  grazing  angles  between  25°  and  65°.  However,  this  model  still 
needs  to  be  verified  with  backscattered  signal  acquired  near  normal  incidence  from  sea 
temperature  around  30°  C. 

On  the  second  task,  we  have  modified  an  existing  stand-alone  L-band  sea  clutter  pro¬ 
gram  written  in  “C”  for  KHOROS  compliance  and  integrated  the  model  with  existing 
RLSTAP  programs  with  minimum  code  modification.  This  exercise  has  allowed  us  to 
examine  the  ability  to  add  functionality  and  to  further  enhance  RLSTAP  modelling  capa¬ 
bilities.  The  existing  sea  clutter  program  in  RLSTAP,  namely  the  “RSS  SigmaO  Curve” 
program  is  inadequate  for  modelling  L-band  VV-polarised  sea  clutter  data.  The  model  is 
not  representative  of  the  measured  NRL-4FR,  JOSS  I,  JOSS  II  and  Darwin  POLSAR  sea 
clutter  data.  The  addition  of  the  MC  sea  clutter  glyph  allows  more  accurate  modelling  of 
the  microwave  returns  from  the  sea  surface  for  radar  operating  in  the  L-band  frequencies. 


Acknowledgements 

The  author  would  like  to  acknowledge  the  contributions  of  Drs.  Nick  Stacy,  and  Ed 
Kruzins  (and  team)  for  providing  the  POLSAR  Darwin  images  and  ground  truth  data. 


25 


DSTO-TR-1206 


(a)  MC  Model  (wind  speed=1.3  m/s) 
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(b)  RSS  SigmaO  Curve  (Sea  State  =  1) 

Figure  5. 4'-  Comparison  of  MC  and  RSS  SigmaO  Curve  sea  clutter  modelling 
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(b)  RSS  SigmaO  Curve  (Sea  State  =  2) 

Figure  5.5:  Comparison  of  MC  and  RSS  SigmaO  Curve  sea  clutter  modelling 
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Appendix  A:  KHOROS  and  Non-KHOROS 

Source  Files 

A.l  Main  Program  for  KHOROS  Generated  MC.c  Source 
File 

/* 

*  Khoros:  $Id$ 

*/ 

#if  ! defined ( lint)  kk  ! def ined ( CODECENTER ) 

static  char  rcsid[]  =  "Khoros:  $Id$"; 

#endif 

/ 

*  Copyright  (C)  1993  -  1997,  Khoral  Research,  Inc.,  ("KRI") . 

*  All  rights  reserved.  See  $BOOTSTRAP/repos/license/License  or  run  klicense. 

*/ 

/*  »»»»»»»»>»»»»»»  <<<<<<<<<<<<<<<<<<<<<<<<<< 

»» 

>»>  Main  program  for  MC 

»» 

»»  Private: 

»»  main 

»» 

»»  Static : 

»»  Public: 

»» 

»»»»»»»»»»>»»»»  <<<<<<<<<<<<<<<<<<<<<<<<<<  */ 


#include  "MC.h" 

clui_inf o_struct  *clui_info  =  NULL; 

void  SeaClutter (void) ;  /*  #######  (a)  INSERT  MC  FUNCTION  PROTOTYPE  ############  */ 

/* - 

I 

I  Routine  Name:  main()  -  Generates  L-band  Sea  Clutter  SigmaO  vs  Grazing  Angles  LUT 

I 

I  Purpose:  main  program  for  MC 

I 

I  Input : 

I  float  clui_info->f_float ;  {float} 

I  int  clui_info->f_flag;  {TRUE  if  -f  specified} 

I 

I  int  clui_inf o->polarisation_toggle ;  {Polarisation  of  co-polarised  transmit  radar  signal  (W  or  HH)} 
I  int  clui_inf o->polarisation_f lag ;  {TRUE  if  -polarisation  specified} 

I 

I  int  clui_inf o->WindDir_toggle ;  {Integer  toggle  selection} 

I  int  clui_info->WindDir_flag;  {TRUE  if  -WindDir  specified} 

I  float  clui_inf o->WindSpeed_float ;  {Measured  wind  speed  above  sea  surface  in  m/s} 
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int  clui_info->WindSpeed_flag;  {TRUE  if  -WindSpeed  specified} 

char  *clui_inf o->OutputData_f ile ;  {SigmaO  versus  grazing  angles} 
int  clui_inf o->OutputData_f lag ;  {TRUE  if  -OutputData  specified} 

Output : 

Returns : 

Written  By:  Poh  Lian  Choong 
Date : 

Modifications : 

- */ 


int  main( 
int  argc , 
char  **argv) 

{ 

/*  -main_variable_list  */ 

/*  -main_variable_list_end  */ 

khoros_init (argc ,  argv,  "MC",  PRODUCT_RELEASE_DATE , 

PRODUCT _RELEASE_NAME,  PRODUCT_RELEASE_VERS ION, 

PRODUCT _RELEASE_MA JOR ,  PRODUCT_RELEASE_MINOR , 
"$MC/objects/kroutine/MC") ; 
kexit_handler(MC_free_args,  NULL); 

/*  -main_get_args_call  */ 

kclui_init ( "MC" ,  "MC" ,  KGENJCROUTINE,  &clui_uis_spec, 

MC_usage_additions ,  MC_get_args, 

MC_free_cirgs) ; 

/*  -main_get_sirgs_call_end  */ 

/*  -main_bef ore_lib_call  */ 

/*  -main_bef ore_lib_call_end  */ 

/*  -main_library_call  */ 

SeaClutterO ;  /*  ##  (b)  INSERT  CALL  TO  L-band  SEA  CLUTTER  PROGRAM/FUNCTION  ###*/ 
/*  -main_library_call_end  */ 

/*  -main_after_lib_call  */ 

/*  -main_af ter_lib_call_end  */ 


kexit (KEXIT_SUCCESS) ; 

} 


/* 


Routine  Name: 
Purpose : 
Input : 
Output : 


MC_usage .additions 

Prints  usage  additions  in  MC.usage  routine 

None 

None 
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I  Written  By:  ghostwriter  -oname  MC 
I  Date :  July  1 ,  1997 

I  Modifications: 

I 

- */ 

void  MC_usage_additions(void) 

{ 

kfprintf  (kstderr ,  "\tGenerates  L-band  Sea  Clutter  SigmaO  vs  Grazing  Angles  LUT\n")  ; 

/*  -usage_additions  */ 

/♦  -usage_additions_end  ♦/ 


} 

/♦ 


Routine  Name : 
Purpose : 
Input : 


MC_free_args 

Frees  CLUI  struct  allocated  in  MC_get_args () 
None 


Output : 
Written  By: 
Date: 
Modifications : 


None 

ghostwriter  -oname  MC 
July  1,  1997 


- */ 

/♦  ARGSUSED  ♦  / 
void 

MC_f ree_args ( 

kexit_status  status, 
kaddr  client_data) 


/♦do  the  wild  and  free  thing  ♦/ 
if  (clui_info  !=  NULL) 

{ 

kfree_and_NULL(clui_info->OutputData_f ile) ; 
kfree_and_NULL(clui_info) ; 

> 

/♦  -free_handler_additions  */ 

/*  -free_handler_additions_end  ♦/ 

} 
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A. 2  KHOROS  Generated  Header  File  -  MC.h 

/* 

*  Khoros:  $Id$ 

*/ 


/* 

*  Copyright  (C)  1993  -  1997,  Khoral  Research,  Inc.,  ("KRI") . 

*  All  rights  reserved.  See  $BOOTSTRAP/repos/license/License  or  run  klicense . 
*/ 

/*  »>»»»»»»»»»»»»»  <<<<<<<<<<<<<<<<<<<<<<<<<< 

»» 

»>>  Purpose:  Include  file  for  MC 

»» 

»»  Written  By: 

»» 

»>>  Modifications: 

»» 

»>>  Date:  July  1,  1997 

»» 

»»»»»»»»»»»»»»>  <<<<<<<<<<<<<<<<<<<<<<<<<<  */ 


#ifndef  _MC_h_ 

#define  _MC_h_ 

* 

*/ 


/* - 

|  #includes 


#include  <mc.h> 

/*  -include_includes  */ 

/*  -include_includes_end  */ 


/* 


#def ines 


* 

*/ 


/*  -include_def ines  */ 

/*  -include_def ines_end  */ 

/* - * 

|  typedefs 

- */ 


typedef  struct  _clui_inf o_struct  { 

/* 

*  float  (required  float) 

*/ 

float  f_float;  /*  float  FLOAT  */ 
int  f_flag;  /*  float  FLAG  */ 

/* 

*  Polarisation  of  co-polarised  transmit  radar  signal  (W  or  HH)  (Required  integer  toggle) 
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*  1  (1) 

*  2  (2), 

*/ 

int  polarisation_toggle;  /*  Polarisation  of  co-polarised  transmit  radar  signal  (W  or  HH)  INT  TOGGLE  */ 
int  polarisation_flag;  /*  Polarisation  of  co-polarised  transmit  radar  signal  (W  or  HH)  FLAG  */ 

/* 

*  Integer  toggle  selection  (Required  integer  toggle) 

*  1  (0) 

*  2  (90) 

*  3  (180), 

*/ 

int  WindDir .toggle;  /*  Integer  toggle  selection  INT  TOGGLE  */ 
int  WindDir.f lag;  /*  Integer  toggle  selection  FLAG  */ 

/* 

*  Measured  wind  speed  above  sea  surface  in  m/s  (required  float) 

*/ 

float  WindSpeed.f loat ;  /*  Measured  wind  speed  above  sea  surface  in  m/s  FLOAT  */ 
int  WindSpeed.f lag;  /*  Measured  wind  speed  above  sea  surface  in  m/s  FLAG  */ 

/* 

*  SigmaO  versus  grazing  singles  (required  outfile) 

*/ 

char  *0utputData_f ile ;  /*  SigmaO  versus  grazing  angles  FILENAME  */ 
int  OutputData.flag;  /*  SigmaO  versus  grazing  angles  FLAG  */ 

}  clui.info.struct; 

/*  -include.typedef s  */ 

/*  -include.typedef s.end  */ 


/* - * 

I  global  variable  declarations 

- 


extern  clui.info.struct  *clui_info; 
extern  uis.struct  clui.uis.spec; 

/*  -include.variables  */ 

/*  -include.variables.end  */ 


/* - * 

I  macros 

- +/ 

/*  -include.macros  */ 

/*  -include.macros.end  */ 

/* - * 

I  routine  definitions 

- 


int  main  PROTO ((int,  char  **)); 

void  MC.get.args  PROTO ( (kaddr)) ; 

void  MC.usage.additions  PROTO ( (void) ) ; 

void  MC.free.args  PROTO ((kexit.status,  kaddr)); 
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/*  -include_routines  */ 

/*  -include_routines_end  */ 

#endif 
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A. 3  Main  MC  Sea  Clutter  Program  -  “SeaClutter.c” 

#include  <stdio.h> 

#include  <stdlib.h> 

#include  <math.h> 

#include  "MC.li"  /*  #########  (c)  INSERT  KHOROS  HEADER  FILE  ###########*/ 

#ifndef  _FCOMPLEX_DECLARE_T_ 

typedef  struct  FCOMPLEX  {float  r,i;>  fcomplex; 

#def ine  _FCOMPLEX_DECLARE_T_ 

#endif  /*  _FCOMPLEX_DECLARE_T_  */ 

static  float  sqrarg; 

#def ine  SQR(a)  ((sqrarg=(a))  — =  0.0  ?  0.0  :  sqrarg* sqrarg) 
static  float  maxargl  ,maxarg2 ; 

#def ine  FMAX(a.b)  (maxargl=(a) ,maxarg2=(b) , (maxargl)  >  (maxarg2)  ?\ 

(maxargl)  :  (maxarg2)) 

#def ine  FREERETURN  {free_vector (fvec, 1 ,n) ;free_vector (xold, 1 ,n) ; \ 

free_vector(p,l,n) ;free_vector (g,l ,n) ;free_matrix(f jac, 1 ,n, 1 ,n) ;\ 
f ree_ivector (indx ,  1  ,n) jreturn; > 


int  nn; 

float  *fvec,w,z; 

#def ine  NR_END  1 
#def ine  FREE_ARG  char* 

#def ine  HORIZONTAL  1 
#def ine  VERTICAL  2 
#def ine  UPWIND  0 
#def ine  DOWNWIND  180 
#def ine  CROSSWIND  90 
#define  Z  1250 

#define  Km2  13.177e4;  /*  sea  surface  tension  and  density  */ 

#define  El  72  /*  L-band  sea  refractivity  real  value  */ 

#define  E2  -59  /*  L-band  sea  refractivity  complex  value  */ 

#define  NAngle  1000 
#define  psiMin  0.1 
#define  psiMax  89.9 
#define  psil  0.1 
#def ine  SPECULAR  60; 

#def ine  pi  3.14159265 
#def ine  MAXITS  200 
#define  T0LF  1.0e-4 
#def ine  T0LMIN  1.0e-6 
#def ine  T0LX  1.0e-7 
#def ine  STPMX  100.0 
#define  N  1 
#define  ALF  1.0e-4 
#def ine  T0LX  1.0e-7 
#def ine  TINY  1.0e-20 
#define  EPS  1.0e-4 

Main  Files - */ 

void  SeaClutter (void) ; 

void  (*nrfuncv) (int  ,  float  *,  float  *) ; 

void  funcv(int,  float  *, float  *) ; 

void  nevrt (float  *,  int  ,  int  *,void  (*vecfunc)  (int,  float  [] ,  float  [))); 
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/* - 

f complex 

f complex 

f complex 

fcomplex 

f complex 

float 

fcomplex 

void 

float 

int 

fcomplex 

float 

void 

void 

void 

void 


-  Utilities  Files 

Complex (float, float) ; 

Cadd(f complex,  fcomplex); 

Csub(fcomplex,  fcomplex); 

Cmul (fcomplex,  fcomplex); 

Cdiv (fcomplex,  fcomplex) ; 

Cabs (fcomplex) ; 

Csqrt (fcomplex) ; 
nrerror(char  *) ; 

♦vector (long, long) ; 

♦ivector (long , long) 

*cvector(long,long) ; 

**matrix(long, long, long, long) ; 
free_vector (float  *, long, long)  ; 
free_ivector (int  ♦, long, long) ; 
free_cvector (fcomplex  *, long, long) ; 
free_matrix(float  long, long, long, long) ; 


*/ 


void  SeaClutter (void)  /*  #############  (d)  CHANGE  mainO  to  a  SUBROUTINE  ############  ♦/ 
Calculate  the  L-band  backscattering  sea  slutter  returns  using  the  MC  Sea  Clutter  model 
(refer  to  report  DST0-TR-0945) 


Copyright:  Poh  Lian  ChoongQDSTO  Salisbury,  Australia,  27  February  2001 
Version  1.0  05  March  2001  modified  for  RLSTAP  compliance 

- •/ 

void  SeaClutter (void) 

{ 

float  *psi,*psiS,*A,k,*K,*K2,v,V,*U,*f ,p,Wl , *W2 , *sigmaB, *sigmaS, *sigma,R2, s2 ; 

float  tmp,tmpl ,tmp2,tmp3; 

float  u.freq; 

int  pol ,  1 ; 

int  i, index, check; 

fcomplex  Er,*R,*P,*Q,S,Ctmp,a,b; 

FILE  *outputfile; 

kfile  *output_f ile ;  /*  ####  (e)  CHANGE  FILE  DECLARATION  to  KH0R0S  FILE  DECLARATION  ####  */ 


/*  ###  (g)  REASSIGN  KH0R0S  CLUI  VARIABLES  to  PROGRAM 


freq 

pol 

u 

1 


=  clui_inf o->f _f loat ; 

=  clui_inf o->polarisation_toggle ; 
=  clui_inf o->WindSpeed_f loat ; 

=  clui_inf o->WindDir_toggle ; 


VARIABLES  ###########  ♦/ 


psi 

psiS 

R 

P 

Q 

A 

K 

K2 

W2 


sigmaB 

sigmaS 

sigma 

U 


=  vector (1 ,NAngle) ; 

=  vector (1 ,NAngle) ; 

=  cvector (1 .NAngle) ; 
=  cvector (1, NAngle) ; 
=  c vect or (1, NAngle) ; 
=  vector (1 .NAngle) ; 

=  vector (1 .NAngle) ; 

=  vector (1 .NAngle) ; 

=  vector (1 .NAngle) ; 

=  vector (1 .NAngle) ; 

=  vector (1 .NAngle) ; 

=  vector (1 .NAngle) ; 
=  vector(l ,2) ; 
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f  =  vector (1,2) ; 

for  (i=l;i<(psiMax*10+l) ;i++)  *(psi+i)=psil*i*pi/180; 

/*  Determine  L-band  sea  reflection  coefficient  */ 

Er  =  Complex(El,E2) ; 

for  (i=l;i<(psiMax*10+l) ;i++){ 

Ctmp  =  Complex(cos(*(psi+i))*cos(*(psi+i)) ,0) ; 

*(R+i)  =  Csqrt(Csub(Er,Ctmp)) ; 

*(P+i)  =  Complex(sin(*(psi+i)) ,0) ; 

*(Q+i)  =  Cmul (Complex (cos (*(psi+i)) ,0) ,Complex(cos(*(psi+i)) ,0)) ;  /*  Q~2  */ 

> 

S  =Csub(Er,Complex(l,0)) ; 

switch(pol){ 
case  HORIZONTAL: 

for  (i=l;i<(psiMax*10+l) ;i++H 
Ctmp  =  Cadd(*(P+i) ,*(R+i)) ; 
b  =  Cmul (Ctmp, Ctmp) ; 

*(A+i)  =  SQR(Cabs (Cdiv(S,b)) ) ; 

> 

break;  ’ 
case  VERTICAL: 

for  (i=l;i<(psiMax*10+l) ;i++){ 

Ctmp  =  Cmul(Er,Cadd(*(Q+i) .Complex (1,0))) ; 
a  =  Cmul(S,Csub(Ctmp,*(Q+i))) ; 

Ctmp  =  Cadd(Cmul(Er,*(P+i)) ,*(R+i)) ; 
b  =  Cmul (Ctmp, Ctmp) ; 

*(A+i)  =  SQR(Cabs(Cdiv(a,b))) ; 

} 

brezik; 

> 


free_cvector(Q,l,NAngle) ; 
free_cvector(P, 1 .NAngle) ; 
free_cvector(R,l,NAngle) ; 

/*  -  Determine  sea  surface  2 

k  =  (2*pi*freq)/3e8; 

for  (i=l;i<(psiHcix*10+l)  ;i++){ 

*(K+i)  =  2*k*cos(*(psi+i)) ; 

* (K2+i)=  SQR(*(K+i))/Km2; 

> 

1  =  1* (pi/180) ; 


•D  directional  sea  spectrum - */ 

/*  radar  wave  number,  k=2*pi/lambda  */ 

/*  Braggs  wave  number  (rad/m)  */ 

/*  Squared  ocean  wave  number  */ 

/*  convert  wind  direction  to  radian  */ 


if  (u  >=  5)  v  =  (0. 003+0. 84e-3*u)/(0. 005+0. 78e-3*u); 

else  v  =  (0.003+1 .92e-3*u)/ (3. 16e-3*u) ;  /*  ratio  of  slope  variance  (clear)  */ 
V  =  2*(l-v)/(l+v) ; 

U[l]=l .0;  z  =  Z,  w=u; 

newt(U,N,A:check,funcv)  ;  /*  Determine  friction  velocity  */ 

funcv(N,U,f ) ; 
p  =  5.0  -  logl0(U[l] )  ; 
f ree_vector (U, 1 ,2) ; 
f ree_vector (f , 1 , 2) ; 
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for  (i=l ; i<(psiMax*10+l) ; i++){ 
tmp2=  *(K+i);  tmp3  =  *(K2+i); 
tmp  =  (1 . 0+3 .0*tmp3) *pow(9. 81 , (l-p)/2) ; 
tmpl=  pow(tmp2,0.922)*pow(tmp2*(l+tmp3) , (l+p)/2) ; 
if  (u<5)  W1  =  1 .750*pow(2*pi ,p-l) *(tmp/tmpl) ; 
else  W1  =  0.875*pow(2*pi,p-l)*(tmp/tmpl) ; 

*(W2+i)  =  (l/(2*pi))*(l+V*cos(2*l))*Wl; 

> 

/*  -  Determine  Bragg’s  component  -  */ 

for  (i=l ; i< (psiMax*10+l) ; i++) { 
tmp  =  SQR(k)*pow(sin(*(psi+i)) ,4) ; 

*(sigmaB+i)  =  10*logl0(tmp* (*(A+i) ) * (* (W2+i) ) ) ; 

* (sigma+i)  =  *(sigmaB+i) ; 

} 

/*  -  Determine  specular  components  -  */ 

a  =  Cmul (Complex (0 . 52,0) ,Csub(Er , Complex (1 ,0) ) ) ; 

b  =  Cmul(Cadd(Csqrt(Er) ,Complex(l , 0) ) ,Cadd (Csqrt (Er) ,Complex(l,0))) ; 

R2  =  SQR(Cabs(Cdiv(a,b))); 

if  (u>=5)  s2  =  1.56e-3*u; 
else  s2  =  5.12e-3*u; 

for (i=601 ; i<(psiMcLX*10+l)  ;i++)  { 

*(psiS+i)=(90. 0  -  i*psil)*pi/180; 
tmpl  =  (R2/s2)*pow(l/cos(*(psiS+i)) ,4) ; 
tmp2  =  exp(-SQR(tan(*(psiS+i)))/s2) ; 

* (sigmaS+i)=10*logl0(tmpl*tmp2) ; 

> 

/*  -  Determine  the  intersection  of  Bragg  and  specular  components  -  */ 

for(i=601 ; i<(psiMax*10+l) ; i++)  { 
if  (*(sigmaS+i)  >=* (sigmaB+i) )  { 
index  =  i ; 
break; 

> 

> 

for (i=(index+l) ; i< (psiMax*10+l) ;i++)  *(sigma+i)=*(sigmaS+i) ; 

/*  ###########  (h)  Change  "C"  FILE  1/0  to  KH0R0S  FILE  1/0  ##############  */ 
output_f ile=kf open(clui_inf o->0utputData_f ile, "w") ; 
for (i=l ; i<=index; i++) 

kf  printf  ( output _f  ile,  "'/,f  7,f  \n"  ,  i*psil ,  *  (sigma+i)  )  ; 
for (i=(index+l) ; i< (psiMax*10+l) ;i++) 

kf  printf  (output_f  ile,  "'/.f  ’/.f  \n"  ,i*psil,  *(  sigma+i) )  ; 
kf close (output_f ile) ; 

free_vector (sigma, l.NAngle) ; 
free_vector(sigmaS,l,NAngle) ; 
free_vector(sigm^, l.NAngle)  ; 
free.vector (W2, l.NAngle) ; 
free_vector (K2, l.NAngle) ; 
f ree_vector (K , 1 , NAngle) ; 
free_vector (A, 1 .NAngle) ; 
free_vector (psiS, l.NAngle) ; 
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free_vector(psi,l,NAngle) ; 

}/*  SeaClutterO  */ 

/* - 

OBJECTIVE  FUNCTIONS 

- */ 

void  funcv(int  n,  float  x[]  .float  f[]) 

< 

float  L; 

L  =  (0.684/x[l])+4.28e-5*SqR(x[l])-0.0443; 
f  Cl]  =  (x[l]/0.4)*log(z/L)*le-2-w; 

} 

/* - 

OPTIMISATION  ROUTINES  (from  Numerical  Recipes  in  C) 

- */ 

void  nevt(float  x[]  ,  int  n,  int  *check,void  (*vecfunc) (int,  float  []  ,  float  [])) 

{ 

void  fdjac(int  n,  float  x[],  float  fvec[],  float  **df,void  (*vecfunc) (int,  float  [] ,  float  [])); 
float  fmin(float  x[]); 

void  lnsrch(int  n,  float  xold[]  .float  fold, float  g[], float  p[]  .float  x[], float 
*f  .float  stpmax,  int  *check,  float  (*func)  (float  [])); 
void  lubksb(float  **a,  int  n,  int  *indx,  float  b[]); 
void  ludcmp(float  **a,  int  n,  int  *indx,  float  *d) ; 


int  i.its, j,*indx; 

float  d.den.f .fold, stpmcix, sum, temp, test, **fjac,*g,*p,*xold; 

indx=ivector (1 ,n) ; 
f jac=matrix(l,n,l,n) ; 
g=vector (1 ,n) ; 
p=vector (1 ,n) ; 
xold=vector ( 1 , n) ; 
fvec=vector (l,n) ; 


nrf uncv=vecf unc ; 
f =fmin(x) ; 
test=0.0; 

for  (i=l ; i<=n; i++) 

if  (fabs(fvec[i])  >  test)  test=f  abs  (fvec  [i]  )  ; 
if  (test  <  0.01+T0LF)  { 

*check=0; 

FREERETURN 

> 

for  (sum=0 .0, i=l ; i<=n;i++)  sum  +=  SQR(x[i]); 
stpmax=STPMX*FMAX(sqrt (sum) , (float)n) ; 
for  (its=l;its<=MAXITS;its++)  { 

fdjac (n,x,fvec,f jac.vecfunc) ; 
for  (i=l;i<=n;i++)  { 

for  (sum=0.0,  j=l;  j<=n;  j++)  sum  +=  f  jac [j]  [i] *fvec [j]  ; 
g[i]=sum; 

> 

for  (i=l;i<=n;i++)  xold[i]=x[i]  ; 
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f old=f ; 

for  (i=l ; i<=n; i++)  p[i]  =  -fvec[i]; 
ludcmp(f jac,n,indx,&d) ; 
lubksb(f jac,n,indx,p) ; 

lnsrch(n,xold,f old, g,p,x,fcf .stpmax, check, f min) ; 
test=0.0; 

for  (i=l ; i<=n;i++) 

if  (fabs(fvec[i])  >  test)  test=f  abs  (f  vec  [i]  )  ; 
if  (test  <  TOLF)  { 

*check=0 ; 

FREERETURN 

} 

if  (*check)  ■( 

test=0.0; 

den=FMAX(f ,0 . 5*n) ; 
for  (i=l ; i<=n; i++)  { 

temp=f abs (g [i] )*FMAX(fabs(x[i] ) , 1 .0)/den; 
if  (temp  >  test)  test=temp; 

> 

*check=(test  <  TOLMIN  ?  1  :  0) ; 

FREERETURN 

> 

test=0.0; 

for  (i=l;i<=n;i++)  { 

temp=(f  abs  (x[i]  -xold[i]  ))/FMAX(f  abs(x  [i]  )  ,1 .0)  ; 
if  (temp  >  test)  test=temp; 

} 

if  (test  <  TOLX)  FREERETURN 

} 

nrerrorC'MAXITS  exceeded  in  newt"); 

} 

/* - 

OPTIMIZATION  UTILITIES  (from  Numerical  Recipes  in  C) 


*/ 


void  lnsrch(int  n,  float  xold[]  ,  float  fold,  float  g[]  ,  float  p[]  ,  float  x[]  , 

float  *f,  float  stpmax,  int  *check,  float  (*func)  (float  [])) 


-c 


int  i ; 

float  a,alam,alam2, alamin,b,disc , f 2 ,f old2 ,rhsl ,rhs2, slope , sum, temp,  test .tmplam; 
*check=0 ; 

for  (sum=0.0,i=i;i<=n;i++)  sum  +=  p[i]*p[i]; 
sum=sqrt(sum) ; 
if  (sum  >  stpmax) 

for  (i=l;i<=n;i++)  p[i]  *=  stpmax/sum; 

for  (slope=0.0,i=l;i<=n;i++)  slope  +=  g[i]*p[i]; 

test=0.0; 

for  (i=l; i<=n; i++)  { 
temp=fabs(p[i] )/FMAX(fabs(xold  [i]  )  ,  1 .0)  ; 
if  (temp  >  test)  test=temp; 

> 

alamin=TOLX/test ; 
alsim=l  .0; 
for  ( ; ;)  { 

for  (i=l;i<=n;i++)  x[i]=xold[i]+alam*p[i]  ; 

*f=(*func) (x) ; 
if  (alam  <  alamin)  { 
for  (i=l;i<=n;i++)  x[i]=xold[i]  ; 

*check=l; 
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return ; 

}  else  if  (*f  <=  f old+ALF*alam*slope)  return; 
else  { 

if  (alam  ==  1.0) 

tmplam  =  -slope/ (2. 0*(*f-f old-slope) ) ; 
else  { 

rhsl  =  *f-f old-alam*slope ; 
rhs2=f 2-f old2-alam2*slope ; 

a=(rhsl/ (alam*alam)-rhs2/ (alam2*alam2))/(alam-alam2) ; 
b=(-alam2*rhsl/ (alam*alam)+alam*rhs2/ (alam2*alam2) )/ (alam-alam2)  ; 
if  (a  ==  0.0)  tmplam  =  -slope/ (2 . 0*b) ; 
else  { 

disc=b*b-3 . 0*a*slope ; 

if  (disc<0.0)  nrerror ("Roundoff  problem  in  lnsrch."); 
else  tmplam=(-b+sqrt(disc))/(3.0*a) ; 

> 

if  (tmplam>0.5*alam)  tmplam=0. 5*alam; 

> 

> 

alam2=alam; 
f2  =  *f; 
f old2=fold; 

alam=FMAX (tmplam, 0 . l*alam) ; 

> 

}•  /*  InsrcbO  */ 

float  fmin(float  x[]) 

int  i; 
float  sum; 

(♦nrfuncv) (nn.x.fvec) ; 

for  (sum=0 .0, i=l ; i<=nn; i++)  sum  +=  SQR(fvec [i] ) ; 
return  0.5*sum; 

}■  /*  fmin()  */ 

void  lubksb (float  **a,  int  n,  int  *indx,  float  b[]) 

{ 

int  i,ii=0,ip, j ; 
float  sum; 

for  (i=l ;i<=n;i++)  { 
ip=indx[i]  ; 
sum=b[ip]  ; 
b[ip]=b[i]  ; 
if  (ii) 

for  (j=ii;  j<=i-l;  j++)  sum  -=  a[i]  [j]  *b[j]  ; 
else  if  (sum)  ii=i; 
b[i]=sum; 

> 

for  (i=n;i>=l;i — )  { 
sum=b  [i]  ; 

for  (j=i+l;  j<=n;  j++)  sum  -=  a[i]  [j] *b[j]  ; 
b[i]=sum/a[i]  [i]  ; 

> 

}/*  lubksb ()  */ 

void  ludcmp (float  **a,  int  n,  int  *indx,  float  *d) 
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int  i,imax,j,k; 
float  big, dum, sum, temp; 
float  *vv; 

vv=vector(l,n) ; 

*d=l .0; 

for  (i=l;i<=n;i++)  { 
big=0.0; 

for  (j=l; j<=n; j++) 

if  ( (temp=fabs(a[i]  [j]  ))  >  big)  big=temp; 
if  (big  ==  0.0)  nrerror ("Singular  matrix  in  routine  ludcmp"); 
vv[i]=l.  0/big; 

} 

for  (j=l ;  j<=n;  j++)  { 

for  (i=l;i<j ;i++)  { 
sum=a[i]  [j]  ; 

for  (k=l  ;k<i  ;k++)  sum  -=  a [i]  [k]  *a [k]  [j]  ; 
a[i]  [j]  =sum; 

> 

big=0.0; 

for  (i=j;i<=n;i++)  { 
sum=a[i]  [j]  ; 
for  (k=l ;k<j ;k++) 

sum  -=  a[i]  [k]*a[k]  [j]  ; 
a[i]  [j]  =sum; 

if  (  (dum=w[i]  *f  abs  (sum)  )  >=  big)  { 
big=dum; 
imax=i ; 

> 

> 

if  (j  !=  imax)  { 

for  (k=l ;k<=n;k++)  { 

dum=a  [imax]  [k]  ; 
a  [imax]  [k]=a[j]  [k]  ; 
a[j]  [k]  =dum; 

> 

*d  =  -(*d) ; 
vv[imax]=vv[j]  ; 

} 

indx  [j]=imax; 

if  (a[j]  [j]  ==  0.0)  a[j]  [j]  =TINY ; 
if  (j  !=  n)  { 

dum=1.0/(a[j]  [j])  ; 

for  (i=j+l ;  i<=n ;  i++)  a[i]  [j]  *=  dum; 

> 

> 

f ree_vector (vv, 1 ,n) ; 

}/*  ludcmp ()  */ 

void  fdjac(int  n, float  x[], float  fvec[]  .float  **df,void  (*vecfunc)  (int, float  [] 
.float  [])) 

int  i.j; 

float  h,temp,*f; 

f=vector(l,n) ; 
for  (j=l; j<=n; j++)  { 
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temp=x[j]  ; 
h=EPS*f abs(temp) ; 
if  (h  ==  0.0)  h=EPS; 
x[j]=temp+h; 
h=x[j]  -temp; 

(♦vecfunc) (n,x,f ) ; 
x  [j]  =temp; 

for  (i=l;i<=n;i++)  df [i] [j] =(f [i] -f vec [i] )/h; 

> 

f ree_vector (f , 1 ,n) ; 

}/*  fdjac()  */ 

/* - 

UTILITIES  (from  Numerical  Recipes  in  C) 

f complex  Complex (float  re,  float  im) 

/*  convert  float  variables  to  complex  numbers  */ 

{ 

f complex  c; 
c .  r=re ; 
c . i=im; 
return  c; 

} 

fcomplex  Cadd(f complex  a,  fcomplex  b) 

/*  add  complex  variables  */ 

{ 

fcomplex  c; 
c.r=a.r+b.r; 
c.i=a.i+b.i; 
return  c; 

} 

fcomplex  Csub (fcomplex  a,  fcomplex  b) 

/*  subtract  complex  variables  */ 

fcomplex  c; 
c.r=a.r-b.r; 
c.i=a.i-b.i; 
return  c; 

> 

fcomplex  Cmul (fcomplex  a,  fcomplex  b) 

/*  multiply  two  complex  variables  */ 

{ 

fcomplex  c; 
c.r=a.r*b.r-a.i*b.i; 
c . i=a . i*b . r+a . r*b . i ; 
return  c; 

} 

fcomplex  Cdiv(f complex  a,  fcomplex  b) 

/*  divide  two  complex  variables  */ 

{ 

fcomplex  c; 
float  r,den; 

if  (fabs(b.r)  >=  fabs(b.i))  { 
r=b.i/b.r; 
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den=b . r+r*b . i ; 
c .r=(a.r+r*a. i)/den; 
c . i=(a. i-r*a. r)/den; 

}  else  { 

r=b.r/b.i; 
den=b.i+r*b.r; 
c .r=(a.r*r+a. i)/den; 
c.i=(a.i*r-a.r)/den; 

} 

return  c ; 


float  Cabs (f complex  z) 

/*  determine  the  absolute  value  of  a  complex  variable  */ 

{ 

float  x,y,ans,temp; 
x=fabs(z.r) ; 
y=fabs(z.i) ; 
if  (x  ==  0.0) 
ans=y ; 

else  if  (y  ==  0.0) 
ans=x; 

else  if  (x  >  y)  { 
temp=y/x ; 

ans=x*sqrt ( 1 . 0+temp*temp) ; 

>  else  { 

temp=x/y ; 

ans=y*sqrt ( 1 . 0+temp*temp) ; 

> 

return  ans; 

> 

f complex  Csqrt(f complex  z) 

/*  determine  the  square  root  of  a  complex  variable  */ 

{ 

f complex  c; 
float  x,y,w,r; 

if  ((z.r  ==  0.0)  (z.i  ==  0.0))  { 
c ,r=0 .0; 
c.i=0.0; 
return  c; 

y  else  { 

x=f abs (z.r); 
y=fabs(z.i)  ; 
if  (x  >=  y)  { 
r=y/x; 

w=sqrt (x)*sqrt (0 . 5* (1 . 0+sqrt (1 .0+r*r)) ) ; 

>  else  { 

r=x/y ; 

w=sqrt(y)*sqrt (0 . 5* (r+sqrt (1 .0+r*r))) ; 

> 

if  (z.r  >=  0.0)  •{ 
c.r=w; 

c.i=z.i/(2.0*v); 

}  else  { 

c.i=(z.i  >=  0)  ?  w  ;  -w; 
c.r=z.i/(2.0*c.i) ; 

} 
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void  nr err or (char  error_text [] ) 

/*  Standard  error  handler  */ 

{ 

fprintf (stderr, "Run-time  error. . An"); 
fprintf  (stderr,  "’/,s\n"  ,error_text) ; 
fprintf (stderr .now  exiting  to  system. . An") ; 
exit(l) ; 

} 

float  *vector(long  nl,  long  nh) 

/*  allocate  a  float  vector  with  subscript  range  v[nl..nh]  */ 
float  *v; 

v=(float  *)malloc((size_t)  ((nh-nl+l+NR_END)*sizeof (float))) ; 
if  (!v)  nrerrorC'allocation  failure  in  vectorO"); 
return  v-nl+NR_END; 

> 

int  *ivector(long  nl,  long  nh) 

/*  allocate  am  int  vector  with  subscript  range  v[nl..nh]  */ 

{ 

int  *v; 

v=(int  *)malloc( (size_t)  ( (nh-nl+l+NR_END)*sizeof (int))) ; 
if  ( !  v)  nrerrorC'allocation  failure  in  ivectorO"); 
return  v-nl+NR_END; 

> 

fcomplex  *cvector (long  nl,  long  nh) 

/*  allocate  a  complex  vector  with  subscript  range  v[nl..nh]  */ 

{ 

fcomplex  *v; 

v= (fcomplex  *)malloc ( (size_t)  ((nh-nl+l+NR_END)*sizeof (fcomplex))) ; 
if  (!v)  nrerrorC'allocation  failure  in  cvectorO"); 
return  v-nl+NR_END; 

> 

float  **matrix(long  nrl,  long  nrh,  long  ncl,  long  nch) 

/*  allocate  a  float  matrix  with  subscript  range  m [nrl. .nrh] [ncl. .nch]  */ 

{ 

long  i,  nrow=nrh-nrl+l,ncol=nch-ncl+i; 
float  **m; 

/*  allocate  pointers  to  rows  */ 

m=(float  **)  malloc((size_t) ((nrow+NR_END)*sizeof (float*))) ; 
if  (!m)  nrerrorC’allocation  failure  1  in  matrixO"); 
m  +=  NR_END; 
m  -=  nrl ; 

/*  allocate  rows  and  set  pointers  to  them  */ 

m [nrl] = (float  *)  malloc( (size.t) ( (nrow*ncol+NR_END)*sizeof (float) )) ; 
if  (!m[nrl])  nrerrorC'allocation  failure  2  in  matrixO"); 
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m[nrl]  +=  NR_END; 
m[nrl]  -=  ncl; 

f or ( i=nrl+l ; i<=nrh ; i++)  m  [i]  =m  [i-1]  +ncol ; 

/*  return  pointer  to  array  of  pointers  to  rows  */ 
return  m; 

> 

void  free_vector (float  *v,  long  nl,  long  nh) 

/*  free  a  float  vector  allocated  with  vector ()  */ 

free ( (FREE_ARG)  (v+nl-NR_END) ) ; 

} 

void  free_ivector (int  *v,  long  nl,  long  nh) 

/*  free  an  int  vector  allocated  with  ivectorO  */ 

{ 

free ( (FREE_ARG)  (v+nl-NR_END) ) ; 

> 

void  free_cvector(f complex  *v,  long  nl,  long  nh) 

/*  free  a  complex  vector  allocated  with  cvectorO  */ 

i 

free ( (FREE.ARG)  (v+nl-NR_END) ) ; 

} 

void  free_matrix (float  **m,  long  nrl,  long  nrh,  long  ncl,  long  nch) 

/*  free  a  float  matrix  allocated  by  matrix ()  */ 

{ 

free ( (FREE_ARG)  (m[nrl] +ncl-NR_END) ) ; 
free ( (FREE.ARG)  (m+nrl-NR_END) ) ; 

} 

/* - 

undef  variables 

- */ 

#undef  HORIZONTAL 
#undef  VERTICAL 
#undef  UPWIND 
#undef  DOWNWIND 
# undef  CROSSWIND 
#undef  Km2 
# undef  FREQ 
#undef  El 
#undef  E2 
#undef  NAngle 
#undef  psiMin 
#undef  psiMax 
#undef  psil 
#undef  SPECULAR 
#undef  MAXITS 
#undef  TOLF 
#undef  TOLMIN 
#undef  TOLX 
#undef  STPMX 
#undef  FREERETURN 
#undef  ALF 
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frundef  TOLX 
frundef  TINY 
frundef  EPS 
frundef  FREERETURN 
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